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Abstract 

We survey the theory and experimental tests for the propagation of cosmic rays in the Galaxy 
up to energies of 10 15 eV. A guide to the previous reviews and essential literature is given, 
followed by an exposition of basic principles. The basic ideas of cosmic-ray propagation are 
described, and the physical origin of its processes are explained. The various techniques for 
computing the observational consequences of the theory are described and contrasted. These 
include analytical and numerical techniques. We present the comparison of models with data 
including direct and indirect - especially gamma-ray - observations, and indicate what we can 
learn about cosmic-ray propagation. Some particular important topics including electrons and 
antiparticles are chosen for discussion. 
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1 Introduction 

Cosmic rays (hereafter CR) are almost unique in astrophysics in that they can 
be directly sampled, not just observed via electromagnetic radiation. Other ex- 
amples are meteorites and Stardust. CR provide us with a detailed elemental 
and isotopic sample of the current (few million years old) interstellar medium 
not available in any other way. It is this which makes the subject especially rich 
and complementary to other disciplines. 

CR have featured frequently in Annual Reviews: about 15 articles from 1952 
to 1989! Reviews include heavy nuclei (jlj), collective transport effects (j2j), com- 
position (0) and propagation Q). Cox's recent Annual Review 'The Three Phase 
Interstellar Medium' (5) contains much discussion of CR as one essential com- 
ponent of the interstellar medium (ISM), but inevitably no mention of their 
propagation. Two recent AnnRev articles (@; Q) give extensive discussion of the 
relation of CR to turbulence which means we do not try to cover this also. 

A basic reference is the book 'Astrophysics of Cosmic Rays' (@) which expounds 
all the essential concepts, and is an update of the classic 'The Origin of Cosmic 
Rays' (jij) which laid the modern foundations of the subject, with an updated 
presentation i n rtioh . Good books for basic expositions are (11; [~L2l ) and for 



high energies (13). A basic text emphasizing theory is (|14l ). while the book 



'Astrophysics of Galactic Cosmic Rays' (|l5l ) gives an valuable overview of the 
experimental data and theoretical ideas as of 2001. The bi-annual International 
Cosmic Ray Conference proceeding^] are also an essential source of information, 
especially for the latest news on the subject. 

Recently a plethora of reviews have appeared on the subject of CR above 10 15 



eV, e.g. (IB; 13; 2); on interactions (flil). on experiments and astrophysics |20l). 



more on astrophysics, propagation and composition |2l|; [z|; I23T), and a review 
of models (0); we recommend (j25l ) and the very up to date (|26l ; 27). Therefore 
this topic has been excluded here. At the lowest energy end, we note that MeV 
particles are non-thermal (even if not relativistic) and must be mentioned in a 



x the recent conferences are available via the NASA Astrophysics Data System ADS 
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review of CR, especially since they are important sources of heating and ionization 
of the interstellar medium |28l). As one example of their far-ranging influence, 
star formation in molecular clouds may be suppressed by CR produced in SNR 
nearby (jjjl ). 

It is worth distinguishing between two ways of approaching CR propagation: 
either from the particle point-of-view, including the spectrum and interactions, 
or treating the CR as a weightless collisionless relativistic gas with pressure and 
energy and considering it alongside other components of the interstellar medium 
(j30l ; I31I). Both ways of looking at the problem are valid up to a point, but for 
consistency a unified approach would be desirable and to our knowledge has never 
been attempted. The nearest approach to this is ((32I; [33I). Most papers address 
exclusively one or the other aspect. The first approach is required for comparison 
with observations of CR (direct and indirect) while the second is required for the 
ISM: stability, heating etc. (@). 

The major recent advances in the field are the high quality measurements 
of isotopic composition and element spectra, and observations by gamma-ray 
telescopes, both satellite and ground-based. Space does not allow any discussion 
of the observational data here, but the figures give an illustrative overview of what 
is now available from both direct (Figs. 1-13) and indirect (7-ray) measurements 
(Figs. 14, 16). Concerning the origin of CR, we follow Cesarsky's 1980 Annual 
Review (jj) 'we will, for the most part, sidestep this problem.' Hence we omit CR 
sources including composition and acceleration; for supernova remnants as CR 
sources the literature can be traced back from the most recent H.E.S.S. TeV 7-ray 
results 113) ■ We also omit solar modulation, galaxy clusters and extragalactic 
CR. We mostly restrict attention to our own Galaxy, but mention important 
information coming from external galaxies (via synchrotron radiation). 

We first introduce the theoretical background, and then consider the confronta- 
tion of theory with observation. A number of particular topics are selected for 
further discussion. 

2 Cosmic-ray Propagation: Theory 
2.1 Basics and Approaches 

We present the basic concepts of CR propagation, and techniques for relat- 
ing these to observational data. Practically all our knowledge of CR propa- 
gation comes via secondary CR, with additional information from 7-rays and 
synchrotron radiation. It is useful at the outset to point out why secondary nu- 
clei in particular are a good probe of CR propagation: the fact that the primary 
nuclei are measured (at least locally) means that the secondary production func- 
tions can be computed from primary spectra, cross-sections and interstellar gas 
densities with reasonable precision; the secondaries can then be 'propagated' and 
compared with observations. 

Since the realization that CR fill the Galaxy it has been clear that nuclear 
interactions imply that their composition contains information on their propaga- 
tion (|44i). A historical event was the arrival of satellite measurements of isotopic 
Li, Be, B in the 1970's |45h. Since then the subject has expanded enormously 
with models of increasing degrees of sophistication. The simple observation that 
the observed composition of CR is different from that of solar, in that rare solar- 
system nuclei like Boron are abundant in CR, proves the importance of propaga- 
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tion in the interstellar medium. The canonical 'few g cm -2 ' of traversed material 
is one of the widest-known facts of cosmic-ray physics. 

At present we believe that the diffusion model with possible inclusion of con- 
vection provides the most adequate description of CR transport in the Galaxy at 
energies below about 10 17 eV so we begin by presenting this model. 

2.2 Propagation equation 

The CR propagation equation for a particular particle species can be written in 
the general form: 



dip(r,p, t) 
di 



q(r,p,t) + V - (D XX V*P-V4>) 



dp pp dp p 2 dp 



P0 - | (V ■ V)iP 



(1) 

Tf T r 



where ijj(f,p,t) is the CR density per unit of total particle momentum p at posi- 
tion r, ip(p)dp = 4irp 2 f(p)dp in terms of phase-space density f(p), q(r,p) is the 
source term including primary, spallation and decay contributions, D xx is the 
spatial diffusion coefficient, V is the convection velocity, diffusive reacceleration 
is described as diffusion in momentum space and is determined by the coefficient 
D pp , p = dp I dt is the momentum gain or loss rate, Tf is the time scale for loss 
by fragmentation, and r r is the time scale for radioactive decay. 

CR sources are usually assumed to be concentrated near the Galactic disk 
and to have a radial distribution like for example supernova remnants (SNR). A 
source injection spectrum and its isotopic composition are required; composition 
is usually initially based on primordial solar but can be determined iteratively 
from the CR data themselves for later comparison with solar. The spallation part 
of q(r,p,t) depends on all progenitor species and their energy-dependent cross- 
sections, and the gas density n(r); it is generally assumed that the spallation 
products have the same kinetic energy per nucleon as the progenitor. K-electron 
capture and electron stripping can be included via Tf and q. D xx is in general a 
function of (r, (3,p/Z) where (3 = v/c and Z is the charge, and p/Z determines the 
gyroradius in a given magnetic field; D xx may be isotropic, or more realistically 
anisotropic, and may be influenced by the CR themselves (e.g. in wave-damping 
models). D pp is related to D xx by D pp D xx oc p 2 , with the proportionality con- 
stant depending on the theory of stochastic reacceleration (|8|; 1461 ) as described 
in Section 2.5. V is a function of r and depends on the nature of the Galactic 
wind. The term in V • V represents adiabatic momentum gain or loss in the 
non-uniform flow of gas with a frozen-in magnetic field whose inhomogeneities 
scatter the CR. Tf depends on the total spallation cross-section and n(r). n{f) 
can be based on surveys of atomic and molecular gas, but can also incorporate 
small-scale variations such as the region of low gas density surrounding the Sun. 
The presence of interstellar helium at about 10% of hydrogen by number must 
be included; heavier components of the ISM are not important for producing CR 
by spallation. This equation only treats continous momentum-loss; catastrophic 
losses can be included via Tf and q. CR electrons, positrons and antiprotons 
propagation constitute just special cases of this equation, differing only in their 
energy losses and production rates. 

The boundary conditions depend on the model; often tp = is assumed at the 
'halo boundary' where particles escape into intergalactic space, but this obviously 
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just an approximation (since the intergalactic flux is not zero) which can be 
relaxed for models with a physical treatment of the boundary. 

Equation (1) is a time-dependent equation; usually the steady-state solution 
is required, which can be obtained either by setting difj/dt = or following the 
time dependence until a steady state is reached; the latter procedure is much 
easier to implement numerically. The time-dependence of q is neglected unless 
effects of nearby recent sources or the stochastic nature of sources are being 
studied. Starting with the solution for the heaviest primaries and using this to 
compute the spallation source for their products, the complete system can be 
solved including secondaries, tertiaries etc. Then the CR spectra at the solar 
position can be compared with direct observations, including solar modulation if 
required. 

Source abundances are determined iteratively, comparing propagation calcula- 
tions with data; for nuclei with very small source abundances, the source values 
are masked by secondaries and cross-section uncertainties and are therefore hard 
to determine. Webber (j47l ) gives a ranking from 'easy' to 'impossible' for the 
possibility of getting the source abundances using ACE data. A recent review of 
the high-precision abundances from ACE is in |4g ) and for Ulysses in (j49|). For a 
useful summary of the various astrophysical abundances relevant to interpreting 
CR abundances see (j50l ). 



2.3 Diffusion 

The concept of CR diffusion explains why energetic charged particles have highly 
isotropic distributions and why they are retained well in the Galaxy. The Galactic 
magnetic field which tangles the trajectories of particles plays a crucial role in 
this process. Typical values of the diffusion coefficient found from fitting to CR 
data is D xx ~ (3 — 5) x 10 28 cm 2 s _1 at energy ~1 GeV/n and it increases with 
magnetic rigidity as i? ' 3 — -R ' 6 in different versions of the empirical diffusion 
model of CR propagation. 

On the "microscopic level" the diffusion of CR results from particle scattering 
on random MHD waves and discontinuities. The effective "collision integral" 
for charged energetic particles moving in a magnetic field with small random 
fluctuations 5B <C B can be taken from the standard quasi-linear theory of 



plasma turbulence (15ll ). The wave-particle interaction is of resonant character 
so that an energetic particle is predominantly scattered by those irregularities 
of magnetic field which have their projection of the wave vector on the average 
magnetic field direction equal to h\ = ±s/ (r g /x), where fi is the particle pitch 
angle. The integers s = 0,1,2... correspond to cyclotron resonances of different 
orders. The efficiency of scattering depends on the polarization of the waves 
and on their distribution in k-space. The first-order resonance s = 1 is the 
most important for the isotropic and also for the one-dimensional distribution 
of random MHD waves along the average magnetic field. In some cases - for 
calculation of scattering at small and for calculation of perpendicular diffusion 
- the broadening of resonances and magnetic mirroring effects should be taken 
into account. The resulting spatial diffusion is strongly anisotropic locally and 
goes predominantly along the magnetic field lines. However, strong fluctuations of 
magnetic field on large scales L ~ 100 pc, where the strength of the random field 
is several times higher than the average field strength, lead to the isotropization 
of global CR diffusion in the Galaxy. The rigorous treatment of this effect is not 
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trivial, since the field is almost static and the strictly one-dimensional diffusion 
along the magnetic field lines does not lead to non-zero diffusion perpendicular 
to B, see ([52) and the references cited there. 

Following several detailed reviews of the theory of CR diffusion |53|; 0; B 
the diffusion coefficient at r g < L can be roughly estimated as D xx ~ 
(5B res / B)~ 2 vr g /3, where 5B res is the amplitude of random field at the resonant 
wave number k res = l/r g . The spectral energy density of interstellar turbulence 
has a power law form w(k)dk ~ k~ 2+a dk, a = 1/3 over a wide range of wave 
numbers 1/(10 20 cm) < k < 1/(10 8 cm), see (0), and the strength of the random 
field at the main scale is 5B ~ 5 //G. This gives an estimate of the diffusion coef- 
ficient D xx « 2 x W 27 (3R]£ cm 2 s" 1 for all CR particles with magnetic rigidities 
R < 10 8 GV, in a fair agreement with the empirical diffusion model (the version 
with distributed reacceleration). The scaling law D xx ~ i? 1 / 3 is determined by 
the value of the exponent a = 1/3, typical for a Kolmogorov spectrum. The- 
oretically (jolt ) the Kolmogorov type spectrum might refer only to some part of 
the MHD turbulence which includes the (Alfvenic) structures strongly elongated 
along the magnetic-field direction and which are not able to provide the significant 
scattering and required diffusion of cosmic rays. In parallel, the more isotropic 
(fast magnetosonic) part of the turbulence, with a smaller value of random field 
at the main scale and with the exponent a = 1/2 typical for the Kraichnan type 
turbulence spectrum, may exist in the interstellar medium (j56l ). The Kraichnan 
spectrum gives a scaling D xx ~ i? 1 / 2 which is close to the high-energy asymptotic 
form of the diffusion coefficient obtained in the 'plain diffusion' version of the em- 
pirical propagation model. Thus the approach based on kinetic theory gives a 
proper estimate of the diffusion coefficient and predicts a power-law dependence 
of diffusion on magnetic rigidity, but the determination of the actual diffusion 
coefficient has to be done with the help of empirical models of CR propagation 
in the Galaxy. 



2.4 Convection 

While the most frequently considered mode of CR transport is diffusion, the ex- 
istence of galactic winds in many galaxies suggests that convective (or advective) 
transport could be important. Winds are common in galaxies and can be CR 



driven (|57|). CR play a dynamical role in galactic halos (|58|; [59j). Convection 
not only transports CR, it can also produce adiabatic energy losses as the wind 
speed increases away from the disk. Convection was first considered by (j60l ) 



and followed up by (|61l : |62| ; 1631 ; 1641 : l65l ). Both 1-zone and 2-zone models have 
been studied: a 1-zone model has convection and diffusion everywhere, a 2-zone 
model has diffusion alone up to some distance from the plane, and diffusion plus 
convection beyond. 

A recent AnnRev on galactic winds (0) does not mention CR, surprisingl 
In the same volume, Cox's article on the Three Phase Interstellar Medium \ 
includes CR as a basic component. Direct evidence for winds in our own Galaxy 
seems to be confined to the Galactic centre region from X-ray images. However 
Cox (jjj) is not sure there is a wind: 'A Galactic wind may be occurring, but I 
do not believe that it carries off a significant fraction of the supernova power 
from the Solar Neighborhood because it would carry off a similar power in the 
pervading cosmic rays...' 

For one-zone diffusion/convection models a good diagnostic is the energy- 



Strong, Moskalenko, and Ptuskin 



7 



dependence of the secondary-to-primary ratio: a purely convective transport 
would have no energy dependence (apart from the velocity-dependence of the 
reaction rate), contrary to what is observed. If the diffusion rate decreases 
with decreasing energy, any convection will eventually take over and cause the 
secondary-to-primary ratio to flatten at low energy: this is observed but con- 
vection ( proposed (|64l ) to explain just this effect) does not reproduce e.g. B/C 



very well (|67|). Another test is provided by radioactive isotopes which effectively 



constrain the wind speed to <10 km s kpc for a speed increasing linearly 



with distance from the disk (|67l ). A value of ~ 15 km s^ 1 (constant speed wind) 
is required to fit B/C even in the presence of reacceleration according to (j68T ) 
which can be compared to 30 km s _1 in the wind model of (j69l ): the latter value 
implies an energy dependence of the diffusion coefficient which may conflict with 
CR anisotropy. 

([H) studied a self-consistent two-zone model with a wind driven by CR and 
thermal gas in a rotating Galaxy. The CR propagation is entirely diffusive in a 
zone \z\ < 1 kpc, and diffusive-convective outside. CR reaching the convective 
zone do not return, so it acts as a halo boundary with height varying with energy 
and Galactocentric radius. It is possible to explain the energy-dependence of the 
secondary-to-primary ratio with this model, and it is also claimed to be consistent 
with radioactive isotopes. The effect of a Galactic wind on the radial CR gradient 
has been investigated (j7ol ) ; they constructed a self-consistent model with the wind 
driven by CR, and with anisotropic diffusion. The convective velocities involved 
in the outer zone are large (100 km s _1 ) but this model is still consistent with 
radioactive CR nuclei which set a much lower limit (j67l ). since this limit is only 
applicable in the inner zone. Observational support of such models would require 
direct evidence for a Galactic wind in the halo. 

2.5 Reacceleration 

In addition to spatial diffusion, the scattering of CR particles on randomly moving 
MHD waves leads to stochastic acceleration which is described in the transport 
equation as diffusion in momentum space with some diffusion coefficient D pp . One 
can estimate it as D pp = p 2 V^ / (9D XX ) where the Alfven velocity V a is introduced 
as a characteristic velocity of weak disturbances propagating in a magnetic field, 
see (0; 14) for rigorous formulas. 

Distributed acceleration in the entire Galactic volume cannot serve as the main 
mechanism of acceleration of CR at least in the energy range 1 — 100 GeV/n. 
In this case the particles of higher energy would spend longer in the system, 
which would result in an increase of the relative abundance of secondary nuclei 
as energy increases, contrary to observation. This argument does not hold at 
low energies where distributed acceleration may be strong and it may explain the 
existence of peaks in the ratios of secondary to primary nuclei at about 1 GeV/n 
if the distributed acceleration becomes significant at this energy. The process 
of distributed acceleration in the interstellar medium is also referred to as 'reac- 
celeration' to distinguish it from the primary acceleration process which occurs 



in the CR sources. It has been shown ((7l|; 1461 ) that the observed dependence 
of abundance of secondary nuclei on energy can be explained in the model with 
reacceleration if the CR diffusion coefficient varies as a single power law of rigidity 
D xx ~ R a with an exponent a ~ 0.3 over the whole energy range (corresponding 
to particle scattering on MHD turbulence with a Kolmogorov spectrum), and if 
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the Alfven velocity is ~ 30 km s , which is close to its actual value in the 
interstellar medium. 

In addition to stable secondary nuclei, the secondary K-capture isotopes are 
useful for the study of possible reacceleration in the interstellar medium (0). 
The isotopes 37 Ar, 44 Ti, 49 V, 51 Cr and some others decay rapidly by electron 
capture at low energies where energetic ions can have an orbital electron. The 
probability to have an orbital electron depends strongly on energy and because 
of this the abundance of these isotopes and of their decay products are strong 
functions of energy and sensitive to changes of particle energy in the interstellar 
medium. The first measurements of an energy-dependent decay of V and 51 Cr 
in CR |73l ) were used to test the rate of distributed interstellar reacceleration fTl) 
but refinement of nuclear production cross sections is required to draw definite 
conclusions. 

The gain of particle energy in the process of reacceleration is accompanied by 
a correspon ding energy loss of the interstellar MHD turbulence. According to 
calculations (|75l ) the dissipation on CR may significantly influence the Kraichnan 
nonlinear cascade of waves at less 10 13 cm and even terminate the cascade at 
small scales. This results in a self-consistent change of the rigidity dependence 
of diffusion coefficient with a steep rise of D xx to small rigidities. The scheme 
explains the high-energy scaling of diffusion D xx ~ R 5 and offers an explanation 
of the observed energy dependence of primary to secondary ratios. 

As mentioned above, the data on secondary nuclei provide evidence against 
strong reacceleration in the entire Galaxy at 1 — 100 GeV/n. However, the spectra 
of secondaries can be considerably modified due to processes in the source regions, 
with a small total Galactic volume filling factor, for the regions where the high- 
velocity SNR shocks accelerate primary CR. Two effects could be operating there 
and both lead to the production of a component of secondaries with flat energy 



spectra (|76l : 1 771 ). One effect is the production of secondaries in SNR by the 
spallation of primary nuclei which have a flat source energy spectrum close to E~ 2 . 
Another effect is the direct acceleration by strong SNR shocks of background 
secondary nuclei residing in the interstellar medium; again the secondaries acquire 
the flat source energy spectrum. Calculations (|77l ) showed that these effects 
might produce a flat component of secondary nuclei rising above the standard 
steep spectrum of secondaries at energies above about 100 GeV/n. 



2.6 Galactic structure 

Almost all aspects of Galactic structure affect CR propagation, but the most 
important are the gas content for secondary production and the interstellar ra- 
diation field and magnetic field for electron energy losses. The magnetic field 
is clearly also important for diffusion but the precise absolute magnitude and 
large-scale structure are less important (at least for CR below 10 15 eV) than the 
turbulence properties. 

The distribution of atomic hydrogen is reasonably well known from 21-cm 
surveys, but the molecular hydrogen is less well known since it has to use the CO 
molecular tracer and the conversion factor is hard to determine and may depend 
on position in the Galaxy. In fact CR-gas interactions provide one of the best 
methods to determine the molecular hydrogen content of the Galaxy, because 
of its basic simplicity, as we describe in the chapter on gamma rays. For more 
details we refer to another review (j78l ). 



Strong, Moskalenko, and Ptuskin 



9 



The Galactic magnetic field can be determined from pulsar rotation and dis- 
persion measures combined with a model for the distribution of ionized gas. A 
large-scale field of a few \xG aligned with spiral arms exists, but there is no 
general agreement on the details (|79l ). One recent analysis gives a bisymmetric 
model for the large-scale Galactic magnetic field with reversals on arm-interarm 
boundaries (jiol ). Independent estimates of the strength and distribution of the 
field can be made by simultaneous analysis of radio synchrotron, CR and 7-ray 
data, and these confirm a value of a few [xG, increasing towards the inner Galaxy 
(|8ll ). A lot of effort has gone into constructing magnetic field models to study 
propagation at energies >10 15 eV where the Larmor radius is large enough for the 
global topology to be important; this is relevant to CR anisotropy and the search 
for point sources. Since this is excluded from our review we refer the reader to 

Hi Si!). 

The interstellar radiation field (ISRF) comes from stars of all types and is 
processed by absorption and re-emission by interstellar dust; it extends from the 
far-infrared though optical to the UV. Computing the ISRF is difficult, but a 
great deal of new information on the stellar content of the Galax y a nd dust is 
now available to make better models for use in propagation codes (|85l ). 

The local environment around the Sun (|86l ) is also important, for example the 
local bubble can have an effect on radioactive nuclei as described in Sec 3.2. 
Extensive coverage of the local environment including CR is in |87l). 



2.7 Interactions 

This large subject is well covered in the literature. Details of the essential pro- 
cesses with references have been conveniently collected in our series of papers: 



energy losses of nuclei and electrons (|67l ) , bremsstrahlung and synchrotron emis- 



sion (18ll ). inverse-Compton emission includin g a nisotropic scattering (|88l ). pion 
production of 7-rays, electrons and positrons (|39l ). Pion production has recently 
been studied in great detail using modern particle-physics codes (HI; M), the for- 



mer giving spectra harder by 0.05 in the index and 7-ray yields somewhat higher 
at a few GeV than older treatments. New more accurate parameterizations will 
be important for the new generation of CR and 7-ray experiments. A useful 
guide to spallation cross-section measurements and models is the contribution by 
J. Connell in (jijll ). and a summary of recent advances is in (I92I: I93I : I94] ) . Accounts 
of radioactive and K-capture processes are in 



2.8 Weighted Slabs and Leaky Boxes 

As mentioned at the start of this section, at present we believe that the diffusion 
model with possible inclusion of convection provides the most adequate descrip- 
tion of CR transport in the Galaxy at energies below about 10 17 eV. The closely 
related leaky-box and weighted slab formalisms have provided the basis for most 
of the literature interpreting CR data. 

In the leaky-box model, the diffusion and convection terms are approximated 
by the leakage term with some characteristic escape time of CR from the Galaxy. 
The escape time T esc may be a function of particle energy (momentum), charge, 
and mass number if needed, but it does not depend on the spatial coordinates. 
There are two cases when the leaky box equations can be obtained as a correct 
approximation to the diffusion model: 1) the model with fast CR diffusion in the 
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Galaxy and particle reflection at the CR halo boundaries with some probability 
to escape (0), 2) the formulae for CR density in the Galactic disk in the flat halo 
model (zh <C R) with thin source and gas disks (z gas <C z^) which are formally 
equivalent to the leaky-box model formulae in the case when stable nuclei are 
considered (fiol ). The nuclear fragmentation is actually determined not by the 
escape time T esc but rather by the escape length in g cm 2 : x = vpr esc , where p 
is the average gas density of interstellar gas in a galaxy with the volume of the 
cosmic ray halo included. 

The solution of a system of coupled transport equations for all isotopes involved 
in the process of nuclear fragmentation is required for studying CR propaga- 
tion. A powerful method, the weighted-slab technique, which consists of splitting 
the problem into astrophysical and nuclear parts was suggested for this problem 
([99l : 0) before the modern computer epoch. The nuclear fragmentation prob- 
lem is solved in terms of the slab model wherein the CR beam is allowed to 
traverse a thickness x of the interstellar gas and these solutions are integrated 
over all values of x weighted with a distribution function G(x) der i ved from 
an astrophysical propagation model. In its standard realization ( IPC ; lOll ) the 



weighted-slab method breaks down for low energy CR where one has strong en- 
ergy dependence of nuclear cross sections, strong energy losses, and energy de- 
pendent diffusion. Furthermore, if the diffusion coefficient depends on the nu clear 
species the method has rather significant errors. After some modification (fl~02h 
the weighted-slab method becomes rigorous for the important special case of sep- 
arable dependence of the diffusion coefficient on particle energy (or rigidity) and 
position with no convective transport. The modified weighted-slab method was 
applied to a few simple diffusion models in (jH^ : [t3 ) . The weighted-slab method 
can also be applied to the solution of the leaky-box equations. It can easily be 
shown that the leaky-box model has an exponential distribution of path lengths 
G(x) oc exp(— x/X) with the mean grammage equal to the escape length X. 

In a purely empirical approach, one can try to determine the shape of the 
distribution function G(x) which best fits the data on abundances of stable pri- 
mary and secondary nuclei |l|). It has been established that the shape of G{x) 
is close to exponential: G(x) oc exp(—x/X(R,P)), and this justifies the use of 
the leaky -box model in this case. There are several recent calculations of G{x) 

(fioi Hoi Klz!). 

The possible existence of truncation, a deficit at small path lengths (below a 
few g cm -2 at energies near 1 GeV/n), relati ye to an expo nential path- length 
distribution, has been discussed for decades (Q; [lOll; [lol; [Oil). The problem was 
not solved mainly because of cross-sectional uncertainties. In a consistent theory 
of CR diffusion and nuclear fragmentation in the cloudy interstellar medium, the 
truncation occurs natu rally if some fraction of CR sources resides inside dense 
giant molecular clouds (103). 

For radioactive nuclei, the classical approach is to compute the 'surviving frac- 
tion' which is the ratio of the observed abundance to that expected in the case 
of no decay. Often the result is given in the form of an effective mean gas den- 
sity, to be compared with the average density in the Galaxy, but this density 
should not be taken at face value. The surviving fraction can better be related 
to physical parameters ( 1081 ). None of these methods can face the complexities of 
propagation of CR electrons and positrons with their large energy and spatially 
dependent energy losses. 
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2.9 Explicit models 

Finally the mathematical effort required to put the 3-D Galaxy into a 1-D for- 
malism becomes overwhelming, and it seems better to work in physical space 
from the beginning: this approach is intuitively simple and easy to interpret. 
We can call these 'explicit solutions.' The explicit solution appr oach including 
secondaries was pioneered by and applied to newer data by (jl09l : 1(351 ) with 
analytical solutions for 2D diffusion-convection models with a cosmic-ray source 
distribution, which however had many restrictive approximations to make them 
tractable (no energy losses, simple gas model). More recently a semi-empirical 
model whi ch is 2D and includes energy-losses and reacceleration has been devel- 
oped l68l ). This is a closed- form solution expressed as a Green's function to 



be integrated over the sources. It incorporates a radial CR so urce distribution, 
but the gas model is a simple constant density within the disk. ( llll ) give an an- 



alytical solution for the time-dependent case with a generalized gas distribution 
but now without energy losses. (This shows again the problem of handling both 
gas and energy losses simultaneously in analytical schemes). 

A 'myriad sources model' (Ill2l ). which is act ually a Green's function method 
without energy losses, yields similar results to ( 113l ) for the diffusion coefficient 



and halo size. But applying the no-energy loss case to ACE data is n ot re ally 



justified, and some defects in their formulation have bee n po i nted out (jllll ). A 



3D analytical propagation method has been developed ( 114 ; 1151 ) with energy 



loss and reacceleration, going via a PLD, but it cannot handle ionization losses 
properly, see section 2.2 of ( 1151 ). They have no spatial boundaries at all and 



rather simplified (exponential) forms for the gas and other distributions. An 
appr oach adapted to fine-scale spatial and temporal variations has been described 
( lid ). This uses a Green's function without energy losses or detailed gas model 



and hence is limited in its application, but is useful for studying the effect of 
discrete sources. 

The most advanced explicit solution to date is the fully numerical model de- 
scribed in the next section. Even this has limitations in treating some aspects 
(e.g. when particle trajectories become important at high energies) so one might 
ask whether a fully Monte-Carlo approach (as is commonly done for energies 
> 10 15 eV) would not be better in the future, given increasing computing power. 
This would allow effects like field-line diffusion (important for propagation per- 
pendicular to the Galactic plane) to be explicitly included. However it is still 
challenging: a GeV particle diffusing with a mean free path of 1 pc in a Galaxy 
with 4 kpc halo height takes ~(4000/l) 2 ~ 10 7 scatterings to leave the Galaxy, 
which would even now need supercomputers to obtain adequate statistics. Hence 
we expect numerical solution of the propagation equations to remain an impor- 
tant approach for the foreseeable future. 



2.10 GALPROP 

The GALPROP code (j67l ) was created with the following aims: 1. to enable si- 
multaneous predictions of all relevant observations including CR nuclei, electrons 
and positrons, 7-rays and synchrotron radiation, 2. to overcome the limitations of 
analytical and semi-analytical methods, taking advantage of advances in comput- 
ing power, as CR, 7-ray and other data become more accurate, 3. to incorporate 
current information on Galactic structure and source distributions, 4. to pro- 
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vide a publicly-available code as a basis for further expansion. The first point is 
the most important, the idea being that all data relate to the same system, the 
Galaxy, and one cannot for example allow a model which fits secondary/primary 
ratios while not fitting 7-rays or not being compatible with the known interstel- 
lar gas distribution. There are many simultaneous constraints, and to find one 
model satisfying all of them is a challenge, which in fact has not been met up to 
now. Upcoming missions should benefho 

We give a very brief summary of GALPR OP; for details we refer the reader 



propagation equation (1) is solved numerically on a spatial grid, either in 2D with 
cylindrical symmetry in the Galaxy or in full 3D. The boundaries of the model in 
radius and height, and the grid spacing, are user-definable. In addition there is a 
grid in momentum; momentum (not e.g. kinetic energy) is used because it is the 
natural quantity for propagation in equation (1). Parameters for all processes 
in equation (1) can be controlled on input. The distribution of CR sources can 
be freely chosen, typically to represent SNR. Source spectral shape and isotopic 
composition (relative to protons) are input parameters. Interstellar gas distri- 
butions are based on current HI and CO surveys, and the interstellar radiation 
field is based on a detailed calculation. Cross-sections are based on extensive 
compilations and parameterizations |92l). The numerical solution proceeds in 
time until a steady-state is reached; a time-dependent solution is also an option. 
Starting with the heaviest primary nucleus considered (e.g. 64 Ni) the propagation 
solution is used to compute the source term for its spallation products, which are 
then propagated in turn, and so on down to protons, secondary electrons and 
positrons, and antiprotons. In this way secondaries, tertiaries etc. are included. 
(Production of 10 B via the 10 Be-decay channel is important and requires a second 
iteration of this procedure.) GALPROP includes K-capture and electron strip- 
ping processes, where a nucleus with an electron (H-like) is considered a separate 
species because of the difference in the lifetime. Since H-like atoms have only 
one K-shell electron, the K-capture decay half-life has to be increased by a fac- 
tor of 2 compared to the measured half-life value. Primary electrons are treated 
separately. Normalization of protons, helium and electrons to experimental data 
is provided (all other isotopes are determined by the source composition and 
propagation). 7-rays and synchrotron are computed using interstellar gas data 
(for pion-decay and bremsstrahlung) and the ISRF model (for inverse Compton) . 
Spectra of all species on the chosen grid and the 7-ray and synchrotron skymaps 
are output in a standard astronomical format for comparison with data. Re- 
cent extensions to GALPROP include non-linear wave damping |75l ) and a dark 
matter package. 

The computing resources required by GALPROP are moderate by today's 
standards. We remark that while GALPROP has the ambitious goal of being 
'realistic', it is obvious that any such model can only be a crude approximation 
to reality. Some known limitations are: only energies below 10 15 eV (no trajectory 
calculations), uniform source abundances (no superbubble enhancements), only 
scales >10 pc (no clumpy ISM: limited by computer power), B-field treated as 
random for synchrotron (regular component affects structure of radio emission). 

2 GALPROP has been adopted as the standard for diffuse Galactic 7-ray emission for NASA's 
GLAST 7-ray observatory, and is also made use of by the AMS, ACE, HEAT and Pamela 
collaborations. 

3 http : / / galpr op . st anf or d . edu 



to the relevant papers 
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For these cases other techniques may be more appropriate, and they provide a 
goal for future developments of GALPROP. 



2.11 Numerical versus analytical 

The following expresses the authors' opinion on this matter. The analytical 
approaches are claimed to have various advantages as follows: 

1. Physical insight: of course it is true that analytical solutions for simple 
cases are very useful to get insight into the relations between the quantities 
involved, and for rough estimates. In fact the analytical formulae may become so 
complicated that finally no insight is gained. In contrast the numerical models are 
very intuitive since they generate explicitly the CR distribution over the Galaxy 
for all species. 2. Equivalent to full solution of propagation equation: only 
true under restrictive conditions, especially involving energy- losses and spatially 
varying densities. Electrons and positrons are anyway beyond analytical methods 
(for energy- losses on realistic interstellar radiation fields), while these are an 
essential component of the CR. 3. Faster, easier to compute: with todays 
computers the speed issue has become irrelevant, and the implementation of a 
numerical model is not harder than the complicated integrals over Bessel functions 
etc. 

I n sum mary we can do no better than a quotation from a paper of 26 years ago 
(!) ( 118l ): 'It is unclear whether one would wish to go much beyond the general- 
izations discussed above for an analytically soluble diffusion model. The added 
insight from any analytic solution over a purely numerical approach is quickly 
cancelled by the growing complexity of the formulae. With rapidly developing 
computational capabilities, one could profitably employ numerical solutions...' 
We remark also that for CR air-shower calculations, analytical methods gave 
way to numerical ones at least 40 years ago. 



2.12 Self-consistent models 

A few attempts at a self-consistent description of CR in the Galaxy have been 
made, including them as a relativistic gas as one component of ISM dynamics. 
This is obviously much harder than the phenomenological models described above 
which treat propagation in a prescribed environment. 3D models of the magne- 



tized ISM with a CR-driven wind have been made by (|32j; l33l ) and this is claimed 
to be also consistent with CR secondary/primary ratios. Such a wind has been 
put forward as a possible explanation of the C R g radient problem (|70l ). The 
Park er instability has been re-analysed recently (jl!9l ) using aniso tropic diffusion 
(12(3) and followed by a CR-driven Galactic dynamo model ( 12ll ) which uses an 



extension the Zeus-3D MHD code (|30j) including CR propagation and sources. 
CR propagation in a magnetic field produced by dynamo action of a turbulent 
flow (|3ll ) presents the whole subject from a novel viewpoint. The extension of 
such approaches to include CR spectra, secondaries, 7-rays etc., which would 
provide a complete set of comparisons with observations, would be very desirable 
but has not yet been attempted. Another kind of self-consistency is to include 
the effect of CR on the diffusion coefficient (j75l ) as described in more detail in 
Section 2.5. 
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3 Confrontation of Theory with Data 
3.1 Stable secondary /primary ratios 

The reference ratio is almost always B/C because B is entirely secondary, the 
measurements are better than for other ratios and are available up to 100 GeV. 
Because C,N,0 are the major progenitors of B, th e pr oduction cross sections are 



better known than, e.g. in the case of Be and Li (1221; |94|). 

The usual procedure is to use a leaky-box or weighted-slab formalism with the 
empirical rigidity-dependence X{R) = {(3//3 )X , ((3//3 )(R/R )~ a X for R < 
Rq, R > Rq respectively. The break at Rq is required because B/C is observed 
to decrease to low energies faster than the /3-dependence (which just describes 
the velocity effect on the reaction rate). The source composition depends on the 
form and parameters of X(R), and vice- versa, (since for example B is produced 
by C,N,0 etc) so the procedure is iterative, starting from a solar-like composition. 

A typical parameter set fitting the data (j69h is a = 0.54 , Xq = 11.8 g cm -2 , Rq 
= 4.9 GV/c, with a source spectrum rigidity index -2.35. In principle all other 
secondary /primary ratios should be consistent with the same parameter set. This 
is generally found to b e the case. As a state-of-the-art application of the weighted 
slab technique of we again refer to (6^). This is applicable to stable nuclei 



only but includes energy losses and gains subject to the limitations described 
in Section 2.8. They apply the method to 1-D disk-halo diffusion, convection, 
turbulent diffusion and reacceleration models cast in weighted-slab form. Fig. [8] 
shows B/C and (Sc+Ti+V)/Fe, so-called sub-Fe/Fe, from their paper. Clearly 
the models cannot be distinguished based on these types of data alone, and they 
can all provide an adequate fit; this shows the importance of using other species 
as well as the ones used here. The models can be used to obtain the injection 
spectrum of primaries, and they find an index of 2.3 to 2.4 for C and Fe in the 
energy range 0.5 - 100 TeV, with the propagated spectrum and data shown in 
Fig. [3 

It has been claimed that no break in X(R) is required to fit Voyager 2 outer- 
heliosphere B/C, N/O and sub-Fe/Fe data extende d to 1.5 GeV, plus HEA03 
data, and adopting suitable solar modulation levels ( 1231 ). Voyager 2 provides a 
unique dataset because of the low solar modulation. 

We consider now explicit models in the sense of Section 2.9. Basically the same 
procedure is adopted, with D XX (R) replacing X{R). Again an ad-hoc break in 
D XX (R) is required in the absence of other mechanisms. Because of the unphysical 
nature of such a D XX (R), many attempts to find a better explanation have been 
made, including convection, reacceleration/ wave damping, and local sources. 

Convection implies an energy-independent escape from the Galaxy so that it 
dominates at low energies as the diffusion rate decreases, so it gives a simple low- 
energy asymptotic X(R) oc (3. However this does not resemble the observed B/C 
energy-dependence, being too monotonic (67). Furthermore quite severe limits 
on the convection velocity come from unstable nuclei. 

Reacceleration affects the secondary-primary ratio as described in Section 2.5. 
Many papers have shown how B/C and other ratios can be reproduced with 
reacceleration at a plausible level and no ad- hoc break in the diffusion coefficient. 
An example is shown in Fig. [H An application including recent ACE (Li, Be, B, 
C) data is in (|94l ). Since reacceleration at some level must be present if diffusion 
occurs on moving scatterers (e.g. Alfven waves) this mechanism is favoured but 
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it is not proven. Direct evidence for reacceleration could come from certain 
K-capture nuclei (Section 3.4). Reacceleration requires a smaller value of a, 
typically 0.3 - 0.4 , consistent with Kolmogorov turbulence, which helps solve 
the problems with anisotropy (section 3.5). 

Closely related to reacceleration is wave-damping as described in Section 2.5. 
This can reproduce B/C, protons and antiprotons satisfactorily as shown in 
Fig. [TU1 and also other data (J75l ). The result of this process is a very sharp rise 
of the diffusion coefficient at rigidities less than about 1.5 GV. The Kolmogorov- 
type dependence is not very successful in this scheme, while a Kraichnan-type 
works better with a high rigidity asymptotic D ~ R 5 (first panel in Fig. [I~0j) . 

A quite different set of parameters has been proposed (68): a = 0.7 — 0.9 
and injection index ~ 2.0, based on fitting many species simultaneously, and is 
suggested to produce the B/C low-energy decrease from convection. Such a large 
a would give problems for the anisotropy (Section 3.5). A related analysis ( lid ) 
claims to exclude the Kolmogorov-type spectrum. 

An alternative explanation for the falloff in B/C at low energies invokes weakly 
nonlinear (in contrast to quasi-line ar) transport theory of cosmic rays in turbulent 
Galactic magnetic fields EH; [Hi). 



Anot her, simp ler, explanation of the B/C energy-dependence is the local-source 
model ( 104 ; 1261 ) in which part of the primary CR have an additional local com- 



ponent. Since the secondary flux has to come from the Galaxy-at-large (the local 
secondaries being negligible) a steep local primary source will cause B/C to de- 
crease at low energies. The known existence of the local bubble containing the 
Sun, and its probable origin in a few supernovae in the last few million years, 
make this a plausible possibility, but hard to prove or disprove. It is claimed 
(Ho!) that if B/C is fitted in such a model, the n sub -Fe/Fe is not fitted, but an 



acceptable fit to this and other data is found in ( 1261 ) using a diffusion model for 
the large-scale component. 

3.2 Unstable secondary /primary ratios: 'radioactive clocks' 

The five unstable secondary nuclei which live long enough to be useful probes 
of cosmic-ray p ropagation ar e 14 C, 10 Be, 26 Al, 36 C1 and 54 Mn, with properties 
summarized in Be is the longest-lived and best measured. The 



theory is pre sented in section 2.2. Based on these isotopes and updated cross- 
sections ( 1281 ) find z&=4-6 kpc, consistent with their earlier estimates of 3-7 kpc 
H) and 4-12 kpc ([67]). Fig. [TT] shows a comparison of 10 Be/ 9 Be with models; 



the ISOMAX 10 Be measurements ( 1291 ) up to 2 GeV (and hence longer decay 



lifetime) are consistent with the fit to the other data, although the statistics are 
not very constraining. 

The data are often inter preted in terms of the leaky-box model, but this is 
misleading ( 1081 : 131 : 127 ). For the formu lae and the detailed procedure for 



the leaky-box model interpretation see (|132l ). Luckily the leaky-box model sur 



viving fraction can be converted to a physically meaningful quantities ( 13ll ) 



for a given model; for example in a simple diffusive halo model, the surviv- 
ing fraction determines the diffusion coefficient, which can be combined with 
stable secondary /primary ratios to derive the halo size. Typical results are 
D xx = (3 — 5) x 10 28 cm 2 s -1 (at 3 GV) and = 4 kpc. We can then compare 
the leaky-box model 'escape time' ~ 10 7 yr with the actual time for CR to reach 
the halo boundary after leaving their sources, the latter being typically an order 
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of magnitude larger. The leaky-box model 'gas density' is typically 0.3 cm~ 3 
compared to the actual average density 0.03 cm -3 for a 4 kpc halo height, again 
an order-of-magnitude difference. 

Since radioactive se cond aries only travel a few hundred pc before decaying it is 
sometimes considered (133) that they cannot give information on the propagation 



region of CR; this is somewhat misleading since it is precisely the combination 
of stable and radioactive data which does allow this; the radioactives determine 
the diffusion coefficient, which then allows the size of the full propagation region 
to be determined from the stable secondary/primary ratio (where the CR do 
come from the entire containment region). This of course assumes the diffusion 
coefficient does not vary greatly from the local region to the full volume. 

The effect of the local bubble surro und ing the Sun on the interpretation of 
radioactive nuclei was pointed out by (|127l ). The flux of unstable secondaries is 
reduced if they are not produced in a gas-depleted region around the Sun, since 
they will decay before reaching us, and this could lead to an overestimate of the 
halo size if interpreted in a simple diffusive halo model. This effect would be 
reduced if the diffusion coefficient in the local region were larger than the large- 
scale value. In fact the situation is even more complex. According to (86; 87; 1341 ) 



the Sun left the local bubble about 10 5 years ago after spending several million 
years inside, and we now live in the CLIC (collection of local interstellar clouds) 
with HI density about 0.2 cm -3 and 35 pc extent. This aspect of the problem 
for CR propagation has not yet been addressed. 

3.3 K-capture isotopes and acceleration delay 

Three isotopes produced in explosive nucleosynthesis decay essentially only by 
K-capture: 59 Ni (7.6 x 10 4 y), 57 Co (0.74 yr), 56 Ni (6d). If acceleration occurs 
before decay, the decay will be suppressed since the nuclei are s tripp ed. 56 Ni is 
absent as expected, but the other two nuclei are more interesting. 

J5a) used ACE 



data on these nuclei to show that the delay between synthesis and acceleration is 
long compared with the 59 Ni decay-time, unless significant 59 Co is synthesised in 
supernovae; considering theoretical 59 Co yields they conclude on a delay > 10 5 y. 
This is inconsistent with models in which supernova accelerate their own ejecta, 
but consistent with acceleration of existing interstellar material. The possibility 
of in-flight elec tron attachment complicates the analysis however (135). For more 



discussion see ()133l ). A result from TIGER on Co/Ni at 1-5 GeV/n (136) supports 



the acceleration delay also at higher energies. 

3.4 K-capture isotopes and reacceleration 

Analyses using ACE data for 49 V 51 V, 51 Cr, 52 Cr and 49 Ti and other nuclei 
have been made Q), (fl37l ). (0). Jl37l ) found that while 51 V/ 52 Cr was in better 



agreement with reacceleration models, 49 Ti/ 46+47+48 Ti gave the opposite result. 
(174! ) find that while V/Cr ratios are in slightly better agreeme nt w ith models 
including reacceleration, the ratios involving Ti are inconclusive. ( 1371 ) used ACE 



data on 37 Ar, 44 Ti, 49 V, 51 Cr, 55 Fe and 57 Co, with inconclusi ye res ults. The main 
problem is the accuracy of the fragmentation cross-sect ions Discussion of 

ACE and previous results can be found in section 3 of ( 1331 ). 
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3.5 Anisotropy 

High isotropy is a distinctive quality of Galactic CR observed at the Earth. The 
global leakage of CR from the Galaxy and the contribution of individual sources 
lead to anisotropy but the trajectories of energetic charged particles are highly 
tangled by regular and stochastic interstellar magnetic fields which isotropize 
the CR angular distribution. This makes difficult or even impossible the direct 
association of detected CR particles with their sources, except for the highest 
energy particles. Observations give the amplitude of the first angular harmonic 
of anisotropy at the level of S ~ 10~ 3 in t he en e rgy range 10 12 to 10 14 eV where 
the most reliable data are available, see (fl~38l : UM and Fig. [T2J The angular 



distribution of particles at lower energies is significantly modulated by the so- 
lar wind. The statistics at higher energies are not good enough yet but the 
measurements indicate the anisotropy amplitude at a level of a few percent at 
10 16 — 10 18 eV. The data of the Super-Kamiokande-I detector ( 140l ) allowed ac- 



curate two-dimensional ma ppin g of CR anisotropy at 10 13 eV (see also the Tibet 
Air Shower Array results (jl38l )). After correction for atmospheric effects, the 
deviation from the isotropic event rate is 0.1% with a statistical significance of 
>5<t and direction to maximum excess at roughly a = 75°, 5 = —5°. 

The amplitude of anisotropy in the diffusion approximation is given by the 
following equation which includes pure diffusion and convection terms: 5 = 
— [3DV/ + up (df /dp)]/vf , see H). Here D is the diffusion tensor, and it is as- 
sumed that the magnetic inhomogeneities which scatter CR particles are frozen 
in the background medium, moving with velocity u <C v , which gives rise to the 
convection term (also called the Compton-Getting term). 

The Compton-Getting anisotropy is equal to (7 + 2)u/c for ultra-rerelativistic 
CR with a power law spectrum 1(E) ~ p 2 f(p) ~ -E~ 7 . The motion of the Solar 
System through the local interstellar medium produces the constant term in the 
energy-dependence of the anisotropy ~4 x 10~ 4 , with a maximum intensity in 
the general direction to the Galactic centre region, which does not agree with 
the data at 10 12 — 10 eV which point to an excess intensity from the anticenter 
hemisphere. The convection effect is outweighed by the diffusion anisotropy which 
is due to the non-uniform distribution of CR in the Galaxy. The systematic 
decrease of CR intensity to the periphery of the Galaxy and the CR fluctuations 
produced by nearby SNR are approximately equally important for the formation 
of the local CR gradient. 

Calculations of CR anisotropy ( 14ll ) are illustrated in Fig. [T2] where the effect 



of global leakage from the Galaxy and the overall contribution of known SNR 
with distances up to 1 kpc are shown separately. Two basic versions of the flat- 
halo diffusion model - the plain diffusion model with D ~ £) - 54 and the model 
with reacceleration where D xx ~ E ' 3 - were used in the calculations. The values 
of D xx were taken from (jH^). The diffusion from a few nearby SNR with different 
ages and distances results in a non-monotonic dependence of anisotropy on energy. 
The flux from the Vela SNR probably dominates over other SNR at energies 
below about 6 x 10 13 eV in the plain diffusion model and below about 10 eV for 
the model with reacceleration. The results of these calculations indicate that the 
diffusion model with reacceleration is compatible with the data on CR anisotropy 
within a factor of about 3. On the other hand it seems that the plain diffusion 
model with its relatively strong dependence of diffusion on energy predicts a too 
large anisotropy at E > 10 14 eV. It should be stressed that CR diffusion is here 
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assumed to be isotropic, which is a seriously simplifying assumption but which 
is difficult to avoid because of the complicated and in fact unknown detailed 
structure of the Galactic magnetic field. The presence of a large-scale random 
magnetic field justifies the approximation of isotropic diffusion on scales larger 
than a few hundred parsecs, but the anisotropy is a characteristic which is very 
sensitive to the local surroundings of the Solar System including the direction of 
magnetic field and the value of diffusion tensor, see discussion in (jl). 



3.6 Diffuse Galactic Gamma Rays 

Gamma rays (above about 100 MeV) from the interstellar medium hold great 
promise for CR studies since they originate throughout the Galaxy, not just the 
local region of direct measurements. The complementarity of 7-rays and direct 
measurements can be exploited to learn most about CR origin and propagation. 
Despite this the interpretation has brought surprises in that the 7-ray spectrum 
is not just as would be expected from the directly-observed CR spectra. 7- 
rays are produced in the interstellar medium by interactions of CR protons and 
He (7r°-decay) and electrons (bremsstrahlung) with gas, and electrons with the 
interstellar radiation field via inverse-Compton scattering. For details of the 
processes the reader is referred to djj§ HH; H; Hljjh . Additional astronomical 
material comes into play like the distribution of atomic and molecular gas, and 
the interstellar radiation field. In fact 7-rays provide an important independent 
handle on molecular hydrogen and its relation to its CO molecular tracer, which 
has taken its place beside more traditional determinations. 

Historically observations started with the OSO-III satellite in 1968, followed 
by SAS-2 in 1972, COS-B (1975-1982) and CGRO (1991-2000). Each of these 
experiments represented a significant leap forward with respect to its predeces- 
sor. SAS-2 established the existence of emission from the interstellar medium 
and allowed a first attempt to derive the CR distribution. With COS-B the CR 
distribution could be better derived and it was found not to follow the canonical 
distribution of supernova remnants, which posed a problem for the SNR origin 
of CR. A rather dependable value for the CO-H2 relation was derived. With 
EGRET and COMPTEL on the Compton Gamma Ray Observatory (CGRO) 
the improvement in data quality was sufficient to allow such studies to be per- 
formed in much greater detail. There is now so much relevant information and 
theory that a rather 'realistic' approach is justified and indeed necessary. At this 
point the best approach seems to be explicit modelling of the high-energy Galaxy 
putting in concepts from CR sources and propagation, Galactic structure, etc. 
The idea of a single model to reproduce both CR and gamma-ray (and other) 
data simultaneously arose naturally a nd is the goal of the GALPROP project 
(see Sec 2.10). For a recent review see ( 1431 ) . 



To illustrate the current state of the art, we show spectra and profiles from a 
recent GALPROP model compared with CGRO /EGRET and COMPTEL data. 
This is based on ( 1131 ) and (|l44i ) where full details can be found. The first model is 



based simply on the directly-measured CR spectra together with a radial gradient 
in the CR sources. It is immediately clear that the spectrum is not well predicted, 
being below the EGRET data for energies above 1 GeV. However remembering 
that this is an unfitted prediction, it does show that the basic assumption that 
the 7-rays are produced in CR interactions is correct; the factor 2 differences are 
telling us something about the remaining uncertainties. The extent to which the 
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CR spectra have to be modified to get a good fit including the GeV excess is shown 
in Fig. [13] and the resulting 7-ray model prediction in Fig. [T4l The difference 
between the directly-observed and modified spectra are within plausible limits 
considering solar modulation and spatial fluctuations in the Galaxy on scales 
>100 pc. A detailed justification is lacking however, so this so-called 'optimized' 
model is just an existence proof rather than a conclusive result. Other more 
drastic modifications to the CR spectrum have been proposed (J78l ) as follows: 

(i) a very hard electron injection spectrum, which could be allowed invoking 
large fluctuations due to energy losses and the stochastic nature of supernova 
remnants in space and time - the solar vicinity is no t necessarily a 'typical' place 
for electron measurements to be representative ( 145 ; 8ll ). Ho weve r the variations 
required are even larger than can reasonably be expected (jll3l ). so this model 
seems unlikely. 

(ii) a hard proton spectrum, again invoking spatial variations in the Galaxy 
so the solar vicinity is not typical; this is more difficult than for electrons since 
the proton energy- losses are negligible. It turns out this possibility can be ruled 
out on the basis of antiproton measurements: too many antiprotons wou ld be 
produced by the same protons which generate the gamma rays ( 1461 : 1131 : l78l ). 
A related suggestion invokes the dispersion in the radio spectral indices of SNR, 
which in dicat es a dispersion in electron indices and if assumed to apply to CR 
protons ( 147l ) could produce the GeV excess; this should also be tested against 
antiprotons. 

Perhaps a completely different source of the excess GeV gam ma r ays is present, 
and the possibility of a dark-matter ori gin has been pursued ( 1481 ) but found to 
produces an excess of CR antiprotons (jl49l ). We will not enter this debate here 
since it is not directly related to the problem of CR propagation. It does however 
show how essential is a good understanding of the CR-induced 7-rays in the 
Galaxy for the study of potentially more fundamental physics ( 150l ). 

The angular distribution of 7-rays provides an essential test for any model. The 
problem with the large expected gradient from SNR is critical. The distribution 
of SNR is hard to measure because of selection effects, so this problem could be 
safely ignored in the past, but now the distribution of pulsars can be determined 
with reasonable accuracy and this should trace SNR since supernovae are pulsar 
progenitors. The pulsar gradient is indeed larger than originally deduced from 
7-rays as shown in Fig. [TBI The distribution of SNR in external galaxies shows a 
similar gradient to pulsars ( 1521 ). The problem therefore remains, but a possible 
solution may be found in another uncertainty, the Galactic distribution of molec- 
ular hydrogen. Since we have to rely on the CO tracer molecule, any vari ation 
in the relation of CO to H2 will affect the interpretation of the 7-ray data. (144]) 
noted that there is independent evidence for an increase in the ratio of H2 to CO 
with Galactocentric radius, related to Galactic metallicity gradients, and that 
this can resolve the problem, allowing CR sources to follow the SNR distribu- 
tion as traced by pulsars. Fig. [16] shows longitude and latitude profiles based on 
such a model, showing a satisfactory agreement with EGRET data. However the 
magnitude of the variation in H2 to CO, or its relation to metallicity and other 
effects, is uncertain so the issue will need more study. Future data from GLAST 
will help by giving much finer angular resolution and the possibility of better 
separating the molecular and atomic hydrogen components. Another possible 
solution to the gradient puzzle in terms of a radially dependent Galactic wind 
was proposed by (|70l ) (see Sections 2.4 and 2.12). 
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Note the large contribution inverse-Compton emission to the spectra and pro- 
files (Figs. [U [TBI) ; this is the reason that the predicted emission is in good 
agreement with the EGRET data right up to the highest latitudes (where the 
gas-relation pion-decay emission is small). Therefore 7-rays constrain both pro- 
tons and electrons, the different angular distributions being clearly distinguish- 
able. The high-latitude inverse-Compton emission is independent evidence for 
the existence of a CR halo extending up to several kpc above the plane, de- 
duced from radioactive CR isotopes as explained in Section 3.2. Interestingly, 
secondary electrons and positrons make a significant contri butio n to the lower- 
energy bremsstrahlung and inverse-Compton gamma-rays ( 1131 ); so we can in 



principle observe secondaries from all over the Galaxy, complementary to the 
local direct measurements. 

3.7 Antiprotons and positrons 

The spectrum and origin of antiprotons in CR has been a matter of active de- 
bate since the first reported detections in balloon flights (fl53l : flolh. There is a 



consen sus t hat most of the CR antiprotons observed near the Earth are secon- 
daries ( 155l ) . Due to the kinematics of secondary production, the spectrum of 



antiprotons has a unique shape distinguishing it from other cosmic-ray species. 
It peaks at about 2 GeV, decreasing sharply towards lower energies. In addition 
to secondary antiprotons there are possible sources of primary antiprotons; those 
most often discussed are dark matter particle annihilation and evaporation of 
primordial black holes. 

In recent years, new data with large statistics on both low and high energy 
antiproton fluxes have become available (fl~56l : [l57l : EH; EH; EH; EH; EH) thanks 



mostly to continuous improvements of the BESS instrument and its systematic 
launches every 1-2 years. This allows one to test models of CR propagation 
and heliospheric modulation. Additionally, accurate calculation of the secondary 
antiproton flux predicts the "background" for searches for exotic signals such as 
WIMP annihilation. 

Despite numerous efforts and overall agreement on the secondary nature of 
the majority of CR antiprotons, pub lished estimates of the expected flux differ 
significantly (see, e.g., Fig. 3 in ( 1581 )). Calculation of the secondary antiproton 



flux is a complicated task. The major sources of uncertainties are three- fold: (i) 
incomplete knowledge of cross sections for antiproton production, annihilation, 
and scattering, (ii) parameters and models of particle propagation in the Galaxy, 
and (iii) modulation in the heliosphere. While the interstellar antiproton flux 
is affected only by uncertainties in the cross sections and propagation models, 
the final comparison with experiment can only be made after correcting for the 
solar modulation. Besides, the spectra of CR nucleons have been directly mea- 
sured only inside the heliosphere while we need to know the spectrum outside, in 
interstellar space, to comput e the antiproton production rate correctly. 

It has recently been shown ( 117l ) t hat a ccurate antiproton measurements during 



the last solar minimum 1995-1997 ( 1581 ) are inconsistent with existing propaga- 



tion models at the ~40% level at about 2 GeV, while the stated measurement 
uncertainties in this energy range are now ~20%. The conventional models based 
on local CR measurements, simple energy dependence of the diffusion coefficient, 
and uniform CR source spectra throughout the Galaxy fail to reproduce simul- 
taneously both the secondary to primary nuclei ratio and antiproton flux. 
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A reacceleration model designed to match secondary to primary nuclei ratios 
produces too few antiprotons because the diffusion coefficient is too large. Models 
without reacceleration can reproduce the antiproton flux, but cannot explain the 
low-energy decrease in the secondary-to-primary nuclei ratios. To be consistent 
with both, the introduction of breaks in the diffusion coefficient and the injection 
spectrum is required, which would suggest new phenomena in particle acceleration 
and propagation. A solution in terms of propagation models requires a break in 
the diffusion coefficient at a few GeV which has been interpreted as a change in 
the propagation mode ( 117l ). The latter calculation employs a modern steady- 
state drift model of propagation in the heliosphere to predict the proton and 
antiproton fluxes near the Earth for different modulation levels and magnetic 
polarity. 

If our local environment influences the spectrum of CR, then it is possible to 
solve the pro blem by invoking a fresh "unprocessed" nuclei component at low 
energies EH), 

which may be produced in the Local Bubble. The idea is that 
primary CR like C and O have a local low-energy component, while secondary 
CR like B are produced Galaxy-wide over the confinement time of 10-100 Myr. 
In this way an excess of B, which appears when propagation parameters are tuned 
to match the p data, can be eliminated by an additional local source of C (see 
Section 3.1). The model appears to be able to describe a variety of CR data, but 
at the cost of additional parameters. 

A consistent p flux in reacceleration models can be obtained if there are addi- 
tional sources of protons <20 GeV (jl63l ). This energy is above the p production 
threshold and effectively produces p's at <2 GeV. The intensity and spectral 
shape of this component could be derived by combining restrictions from p's and 
diffuse 7-rays. 

It is clear that accurate measurements of p flux are the key to testing cur- 
rent propagation models. If new measurements confirm the p "excess," current 
propagation and/or modulation models will face a challenge. If not - it will be 
evidence that the reacceleration model is currently the best one to describe the 
data. Here we have to await for th e nex t BESS-Polar flight which should bring 
much larger statistics, and Pamela ( 164j ). c urren tly in orbit. 

Positrons in CR were discovered in 1964 (|165l ). The ratio of positrons to elec- 
trons is approximately 0.1 above ~1 GeV as measured near the Earth. Secondary 
positrons in CR are mostly the decay products of charged pions and kaons (vr + , 
K + ) produced in CR interactions with interstellar gas. The calculation of the 
positron production by CR generally agrees with the data and indicates that 
the majority of CR positrons are secondary (fi~6i HI). Some small fraction of 
positrons may also come from sources ( 167 ), such as pulsar winds (|l68l ) or WIMP 
annihilations ( 1691 ) . 

In interstellar space the secondary positron flux below ~1 GeV is comparable to 
the electron flux which makes pos itron s non-negligible contributors to the diffuse 
7-ray emission in the MeV range ( 1131 ). The contribution of secondary positrons 
and electrons in diffuse Galactic 7-ray emission models essentially improves the 
agreement with the data, making secondary positrons effectively detectable in 
7-rays. 

The accuracy of spacecraft and balloon-borne experiments has substantially in- 
creased during the last decad e. T he d ata taken by various instruments including 
AMS, CAPRICE and HEAT (Il70l: IZll: [ml : [i~72h are in agreement with each other 
within the error bars, which are however still quite large. A recent measurement 
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of CR positron fraction by the HEAT instrument (|42l ) indicates an excess (fea- 
ture) between 5 and 7 GeV at the level of ~3cr over the smooth prediction of a 
propagation model (j39l ). When all the HEAT flights are combined, the data show 
an excess above ~8 GeV, with the most significant point at ~8 GeV (42). The 
Pamela satellite ( 1641 ) currently in orbit, will have a very good positron statistics. 



3.8 Electrons and Synchrotron Radiation 

CR electrons require a separate treatment from nuclei because of their rapid en- 
ergy losses and their link with synchrotron radiation. Their low energy density 
compared with nuclei (~ 1%) is not yet understood; standard SNR shock acceler- 
ation does not predict such a ratio and it is normally a free parameter in models 
e.g. (jl73l ). Direct measurements extend from MeV to TeV ; at low energies solar 
modulation is so large that the interstellar fluxes are really unknown. At TeV 
energies the statistics are very poor but new experiments are in progress. Young 
SNR contain TeV electrons as is shown by their non-thermal X-ray emission, and 
energies up to 100 TeV are possible ( 1741 ) . The rapid energy losses mean that 



contributions from old nearby SNR such as Loop 1, Vela, the Cygnus Loop and 
MonoGem could make an important contribution to the local electron spectrum 



above 100 GeV and dominate above 1 TeV, see Fig. d] (|37l ). The rapid energy 
losses of electrons also produce time- and space-dependent effects described in 
the next section. The essential propagation effects are seen in Fig. [TBI at energies 
below 1 GeV the interstellar spectrum is flatter than the injection spectrum, due 
to Coulomb losses; there is an intermediate energy range where the spectrum 
is similar to that at injection, and at high energies it steepens due to inverse 
Compton and synchrotron losses and energy-dependent diffusion. For propaga- 
tion models we need an injection spectrum which reproduces the observations; 
a typical spectral index is 2 .4, similar to nuclei, which produces the observed 



high-energy slope of 3.3 ( 113l ). The synchrotron spectral index is /3„ = (7 — l)/2 



for an electron index 7; the observed f3 u = 0.6 — 1, increasing with frequency, 
implies 7 = 2.4 — 3 increasing with energy, in satisfactory agreement. Detailed 
geometric al models of the Galactic synchrotron emission have been constructed 
(jl75l : 1761 ). A more physical interpretation requires propagation and magnetic 



field modelling as for example in (I8ll ). Combined with 7-ray data this allows the 
magnetic field to be determined as a function of Galactocentric radius indepen- 
dent of other techniques (fill ). Since that work a great deal of new radio data 
has become available both from new ground-based surveys and from balloons 
and satellites for cosmology (e.g. WMAP). Since Galactic synchrotron is an es- 
sential 'foreground' for CMB studies there is a large overlap of interest between 
these fields. Exploitation of these data for CR studies has only just begun but is 
expected to accelerate with the forthcoming Planck mission. 

Radio continuum observations of ot her galax ies provide a complementary view 



of electrons via s ynchrotro n radiation (177; 1781 ). The classical case is the edge-on 



galaxy NGC891 (jl79l : [i~8(ih which has a non-thermal halo extending to several kpc, 



and this helped to give credence to the idea of a large halo in our Galaxy. Mor e 
recent observations of this and other edge-on galaxies confirm large halos (fl8lh . 
The spectral index variations due to energy losses provide a test of propagation in 
principle, although this has not be en ve ry fruitful up to now. This effect can also 
be measured in starburst galaxies (182)/ Face-on galaxies like M51 show marked 
spiral structure in the continuum (|177l ). but this is probably mainly a magnetic- 
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field effect. It is not possible reliably to separate CR and magnetic field variations, 
although this is often attempted assuming equipartition of energy between CR 
and field, which although having no firmly establi shed phys ical basis nevertheless 
gives field values quite similar to other methods (jl83l : 1781 ). 

The tight radio-continuum/far- infrared correlation for galaxies ( 184 : 1851 : Il8d ) 
is worth mentioning here. It holds both within galaxies (down to 100 pc or less) 
and from galaxy to galaxy, and over several orders of magnitude in luminosity; it 
presu mably contains clues to CR origin and propagation. The 'CR calorimeter' 
(|187l ) is the simplest interpretation. Simply relating both CR production and 
UV (reprocessed to far-infrared) to star formation rates is apparently insufficient 
to explain the close correlation, and this has given rise to several different inter- 
pretation s incl uding hydrostatic regulation, without being very conclusive as yet 
(fl~8i [Hi; [Hi). 



3.9 Time- and space-dependent effects 



CR propagation is usually modelled as a spatially smooth, temporally steady- 
state process. Because of the rapid diffusion and long containment time-scales in 
the Galaxy this is usually a sufficient approximation, but there are cases where it 
breaks down. The rapid energy loss of electrons above 100 GeV and the stochastic 
nature of their sources produces spatial and temporal variations. Supernovae are 
stochastic events and each SNR source of CR lasts only perhaps 10 — 10 5 years, 
so that they leave their imprint on the distribution of electrons. This leads 
to large fluctuations in the CR electron density at high ene rgies , so that the 
electron spectrum measured near the Sun may not b e typical (|ll3l ); a statistical 
calculation of this effect can be found in (145; 1 1881 ) and for local SNR in (j37l ). 
These effects are much smaller for nucleons since there are essentially no energy 
losses, but they can still be important ( 1891); fo r the theory of CR fluctuations 
for a galaxy with random S NR event s see ()14ll ) and Section 3.5. Such effects 
can influence the B/C ratio ( 111 : lid ) mainly through variations in the primary 
spect ra. T he local bubble can also have an effect on the energy-dependence of 
B/C ijlifih . Dispersion in CR injection spectra from SNR may cause the locally 
observed spectrum to deviate from the average (fl~47h . 



4 Summary Points list 

Progress in CR research is rapid at the present time with the new generation 
of detectors for both direct and indirect measurements. But still it is hard to 
pin-down particular theories and even now the origin of the nucleonic component 
is not proven (|34l ). although SNR are the leading candidates. 
The main points we want to make are 

• the importance of considering all relevant data, both direct (particles) and 
indirect (7-ray, synchrotron) measurements 

• increase in computing power has made many of the old approximations for 
interpreting CR data unnecessary 

• new high-quality data will require detailed numerical models 
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5 Future Issues 

We end by listing some of the open questions regarding CR propagation which 
might be answered with observations in the near future. 

• the interpretation of the energy-dependence of the secondary/primary ratio, 
requiring accurate measurements at both, low and high energies 

• the size of the propagation region - existence of an extended halo, requiring 
measurements of radioactive species over a broad energy range 

• the relative roles of diffusion, convection and reacceleration 

• the importance of local sources to the primary CR flux 

• the origin of the GeV excess in 7-rays relative to the prediction based on 
locally observed CR 

• the CR source distribution: is it like SNR or not? 

• are positrons and antiprotons explained as secondaries from primary CR or 
is there a - perhap exotic - excess? 

In addition, on the theoretical side, we mention 

• the relation of CR-dynamical models of the Galaxy to CR propagation 
theory 



6 Related Resources 



GALPROP website: http://galprop.stanford.edu, 



GLAST website: http://glast.stanford.edu 



7 Key terms, definitions and acronyms 

Key terms/definitions: 

Direct measurements: measurements on CR themselves with detectors in bal- 
loons and satellites. Indirect measurements: via electromagnetic radiation (7- 
rays and synchrotron) emitted by CR interacting with interstellar matter and 
radiation fields. Primary CR: CR accelerated at sources (for example in SNR). 
Secondary CR: CR generated by spallation of primary CR on interstellar gas. p: 
particle momentum, R: particle magnetic rigidity = pc/Ze, j3: velocity/speed- 
of-light, Zh- height of CR halo in direction perpendicular to Galactic plane 

Acronyms: 

ACE: Advanced Composition Explorer, B/C: cosmic-ray Boron-to-Carbon ra- 
tio, CMB: cosmic microwave background, CGRO: Compton Gamma Ray Obser- 
vatory, CR: cosmic rays, EGRET: Energetic gamma-ray telescope, GALPROP: 
Galactic Propagation code, GLAST: Gamma-ray Large Area Telescope, H.E.S.S.: 
High Energy Stereoscopic System, ISRF: interstellar radiation field, PLD: path- 
length distribution, SNR: supernova remnant 
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Figure 1: Compilation of spectral data 10 — 10 eV for p. He, Fe, combining 
balloon, satellite and ground-based measurements. From (|18l ) and G. Horandel, 
private communication. 
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Figure 2: Preliminary spectra of p, He from ATIC-2, compared with AMS01, 
BESS-TeV and CAPRICE98. Plot from ATIC collaboration (35). The ATIC-2 
data indicate a slighly harder spectrum for He above 1 TeV. 
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Figure 3: Compilation of spectral data for element groups CNO, NeMgSi, Fe 
Jiifl) from HEAO-3, SANRIKU, CRN, SOKOL, JACEE and RUNJOB (upper) 
and of separate even-Z elements from preliminary ATIC-2, HEAO-3, CRN and 
TRACER (lower). Plot from ATIC collaboration |H). The ATIC-2 data suggest 
a hardening above 10 TeV. 
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Figure 4: Measurements of the electron spectrum, including AMS01, 
CAPRICE94, HEAT and SANRIKU, compared with possible contributions of 
distant sources and local supernova remnants, from (|37l ). 




Figure 5: Measurements of the positron spectrum, including data from MASS91, 
AMS01, CAPRICE94 and HEAT, from |H). Propagation calculations for in- 
terstellar (upper curves) and modulated (lower curves) are shown. Dotted, dot- 
dashed: GALPROP Hi); dashed, solid: (0). 
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Cosmic-ray propagation 




Figure 6: Left: Measurements of the electron and positron spectrum by AMS01 
and HEAT, from (41). Right: e+/(e+ + e~) ratio including HEAT and AMS01 
(42). For the ratio a propagation calculation (j39l ) is shown by the solid line, 
and possible contributions from pulsars (dashed) and dark matter annihiliations 
(dash-dot). 
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Figure 7: Measurements of the antiproton spectrum from MASS91, CAPRICE98 
and BESS, compared to a propagation calculation, from (43). The dashed and 
dotted lines illustrate the sensitivity of the calculated (modulated) antiproton flux 
to the normalization of the diffusion coefficient. LIS marks the local interstellar 
spectrum. 
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Figure 8: B/C and sub-Fe/Fe data compilation compared to four models treated 
by the modified weighted-slab technique, from (j(59l ). 
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Figure 9: Data compilation for spectra of C and Fe compared to four models 
treated by the modified weighted-slab technique, from (|69i ). 



36 



Cosmic-ray propagation 



10 10 10 10 10 10 
Rigidity, GV 





PROTONS 




DRD model 




* = 550 MV " 










/ O AMS I 




-mj V BESS 98 




- C T • CAPRICE 94 




J? □ IMAX 92 




y A LEAP 87 


ft x Sokol, JACEE 





10 10 10 10 10 10 

Kinetic energy, GeV 



0.3 
0.25 
0.2 



B/C ratio 




DRD model 






* = 450 MV - 








w 


/ \\ 




- LIS/ 


▼ Voyager p 

• Ulysses J. 1 

a ACE 

a HEAO-3 

j Chapell.Webber 1 

■ Dwyer 1 978 

V Maehl etal. 1977 


381 IpL 



10 10 10 10 10 10 

Kinetic energy, GeV/nucleon 

ANTIPROTONS 

• BESS 95-97 
□ BESS 98 
■ MASS 91 
A CAPRICE 




10 10 10 10 

Kinetic energy, GeV 



Figure 10: Upper left: diffusion coefficient in different models (|75l ): plain diffu- 
sion (dots), Kolmogorov reacceleration (dashes), and Kraichnan-type reacceler- 
ation with wave-damping (solid). Upper right: B/C, lower left: protons, lower 
right: antiprotons, in the wave-damping model (|75l ). LIS marks the local inter- 
stellar spectrum. 
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Figure 11: Data on energy-dependence of u Be/ Be including ACE, Ulysses, 
Voyager, IMP, ISEE-3 and ISOMAX data, from (fl29h . The solid line is a diffusive 
halo model with 4 kp c scale height using GALPROP (j9^ ) the other lines are 
leaky-box models ( 130l ). 
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Figure 12: The anisotropy of cosmic rays in the reacceleration (red curves) and 
the plain diffusion (blue curves) models. Separately shown are the effects of the 
global leakage from the Galaxy (thick lines) , and the contribution from local SNR 
(thin lines). The collection of data on cosmic ray anisotropy (yellow squares) are 
taken from ( 1421 ) where the references to individual experiments can be found. 




Figure 13: Directly-observed (solid), and modified (dots) proton and electron 
spect ra. The modified spectra are deduced from fits to antiproton and 7-ray data 
( 1131 ). LIS marks local interstellar spectra. Data shown are from AMS01, BESS, 
CAPRICE94, IMAX, LEAP, HEAT, MASS91, SANRIKU. 
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Cosmic-ray propagation 




Figure 14: 7-ray spectrum of inner Galaxy (330° < I < 30°, \b\ < 5°) for model 
based on the directly-observed CR spectra and modified spectr a sho wn in Fig. [TBI 
Red bars: EGRET data, including points above 10 GeV, see ( 15ll ). Green bars: 
COMPTEL. Light blue line: bremsstrahlung, green line: inverse Compton scat- 
tering, red line: 7r°-decay, black line: extragalac tic b ackground, dark blue line: 
total. This is an update of the spectra shown in ( 1131 ). 
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Figure 15: Possible CR source distributions as function of Galactocentric radius 
R: pulsars (solid line), SNR (vertical bars), from 7-rays assuming constant H2- 
to-CO relation (dotted line) (fui ). 
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Figure 16: 7-ray longitude and latitude prof iles, for model with pulsar source 
distribution and Xco varying with radius (|l44l ). Light blue line: bremsstrahlung, 
green line: inverse Compton scattering, red line: 7r°-decay, black line: extragalac- 
tic background, dark blue line: total. 



